Matrix-free GPU implementation of a preconditioned conjugate gradient solver for anisotropic elliptic PDEs
نویسندگان
چکیده
منابع مشابه
Matrix-free GPU implementation of a preconditioned conjugate gradient solver for anisotropic elliptic PDEs
Many problems in geophysical and atmospheric modelling require the fast solution of elliptic partial differential equations (PDEs) in “flat” three dimensional geometries. In particular, an anisotropic elliptic PDE for the pressure correction has to be solved at every time step in the dynamical core of many numerical weather prediction (NWP) models, and equations of a very similar structure aris...
متن کاملMatrix-free multi-GPU Implementation of Elliptic Solvers for strongly anisotropic PDEs
Manycore chip architectures such as Graphics Processing Units (GPUs) have been shown to be highly e cient both in terms of total performance and small power consumption in many areas of Scientific Computing. However, implementing memory-bound applications such as elliptic solvers for large partial di↵erential equations (PDEs) remains a challenge. We consider iterative solvers for strongly aniso...
متن کاملAn Improved Implementation of Preconditioned Conjugate Gradient Method on GPU
An improved implementation of the Preconditioned Conjugate Gradient method on GPU using CUDA (Compute Unified Device Architecture) is proposed. It aims to solving the Poisson equation arising in liquid animation with high efficiency. We consider the features of the linear system obtained from the Poisson equation and propose an optimization method to solve it. First, a novel storage format call...
متن کاملA fast algebraic multigrid preconditioned conjugate gradient solver
This work presents a new approach for selecting the coarse grids allowing a faster algebraic multigrid (AMG) preconditioned conjugate gradient solver. This approach is based on an appropriate choice of the parameter a considering the matrix density during the coarsening process which implies in a significant reduction in the matrix dimension at all AMG levels. 2005 Elsevier Inc. All rights rese...
متن کاملPetascale elliptic solvers for anisotropic PDEs on GPU clusters
Memory bound applications such as solvers for large sparse systems of equations remain a challenge for GPUs. Fast solvers should be based on numerically efficient algorithms and implemented such that global memory access is minimised. To solve systems with up to one trillion (10) unknowns the code has to make efficient use of several million individual processor cores on large GPU clusters. We ...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
ژورنال
عنوان ژورنال: Computing and Visualization in Science
سال: 2013
ISSN: 1432-9360,1433-0369
DOI: 10.1007/s00791-014-0223-x